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Abstract 

We present and constrain a cosmological model which component is 
a pressureless fluid with bulk viscosity as an explanation for the present 
accelerated expansion of the universe. We study the particular model of 
a constant bulk viscosity coefficient (m- The possible values of are 
constrained using the cosmological tests of SNe la Gold 2006 sample, the 
CMB shift parameter R from the three-year WMAP observations, the 
Baryon Acoustic Oscillation (BAO) peak A from the Sloan Digital Sky 
Survey (SDSS) and the Second Law of Thermodynamics (SLT). It was 
found that this model is in agreement with the SLT using only the SNe 
la test. However when the model is submitted to the three cosmological 
tests together (SNe+CMB+BAO) the results are: 1.- the model violates 
the SLT, 2.- predicts a value of Ho ~ 53 km ■ sec~^ ■ Mpc~^ for the Hub- 
ble constant, and 3.- we obtain a bad fit to data with a Xmin ~ 400 
(Xd.o.f. ~ 2.2). These results indicate that this model is ruled out by the 
observations. 

We present a flat cosmological model which component is a pressureless 
fluid made of baryon and dark matter components with a constant bulk viscosity 
coefficient of the form Cm ^constant, which potentially could explain the present 
accelerated expansion of the universe. A bulk viscosity coefficient can produce 
a positive term in the second Friedmann equation that induces an acceleration 
[T1I2] (i.e., a > , with a denoting the scale factor and the dots derivatives with 
respect to the cosmic time). Similar models or analysis have been proposed also 
in Him [5]. 



Model of the bulk viscous matter-dominated universe The energy- 
momentum tensor of a fluid composed by only matter (hereinafter, we call mat- 
ter to baryon and dark matter components together) with bulk viscosity is 
defined by T^^ = PmUfxU„ -\- (g^^, -\- u^u^)P^, where = - 3Cm (a/a) 
(see e.g. [1]). Here Pm is the pressure of the matter fluid. However, as- 
suming pressureless matter: Pm = ^ Pj* = — 3Cin(a/a)- The subscript 
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"m" stands for matter. The four-velocity vector is that of a comoving 
observer that measures the pressure and the density pm of the matter 
fluid. The conservation of energy has the form a{dpm/da) = 3 (3Cm^^ — Pm), 
where H = a/a is the Hubble parameter. This equation in terms of the red- 
shift reads (1 + z) ^dCl^/dz^ + ClW^ Co ~ 3f),„ = 0, where we have defined 
(ini(z) = pm/p^j.;^., and the dimensionlcss coefficient Co = Cm(247rG/iJo)- The 
exact solution to this equation is Clai{z) = (1/9) (Co — 3)(1 -I- z)^" — Co • On 

the other hand, a spatially flat geometry is assumed (fc = 0) for the FRW metric 
cosmology as suggested by WMAP. Thus, the first Friedmann equation for a flat 
universe composed by bulk viscous matter has the form: H'^{z) ~ (87rG/3)p: 

Hi Cl^iz). Therefore H{z) = (iJo/3) [(Co - s) (1 + zf/^ - Co 



Cosmological tests We constrain the possible values of the bulk viscosity Co 
using the cosmological tests of the Gold 2006 SNe la data sample [6] composed by 
182 SNe la, the Cosmic Microwave Background (CMB) shift parameter R from 
the three-year WMAP observations, the Baryon Acoustic Oscillation (BAO) 
peak A from the Sloan Digital Sky Survey (SDSS) and the Second Law of 
Thermodynamics (SLT). 

For the SNe la we define the observational luminosity distance [7J [S] in a 
flat cosmology and with the bulk viscous matter component as d^, = c(l + 
z)Hq^ E{z')~^ dz' , where E{z) = H{z)/Hq. The theoretical distance moduli 
for the i-th supernovae with redshift Zi is fi^izi) — b\ogif^[dL{zi) /WlYyc] + 25 , 
with c the speed of fight. The Xsno statistical function becomes Xsnc(Co, ^^o) = 

r ~ n 2 

Sfc=i /^*(^fci Co, -f^o) — Pfe /<y\, where pfc is the observed distance moduli for 
the fc-th supernovae and af. is the variance of the measurement. 

The CMB shift parameter R is defined as i? = y/f>^ E{z')-^ dz' , 
where ^cmb = 1089 is the redshift of recombination [Hllin]- The observed value 
of the shift parameter R is reported to be i?obs = 1-70 ± 0.03 [TIT. 

From the LSS data, we use the baryon acoustic oscillation (BAO) peak A 

defined as A = y/^^^E{zi)-^'^ bi^ C E{z')-^ dz']^'"^ , where zi = 0.35 [H]. 
The observed value of SDSS-BAO peak is = 0.469 ± 0.017 [H]. 

We construct the joint function as xLai = xIno + Xcmb + XbaO' where 
XsNc was defined above, and Xcmb = [{R - -Robs)/o-fl]^ and Xbao = " 

Aohs)/(TAY- 

The generation of local entropy in the space-time is defined as TV^s^ = 
(Vfj^u'^ = 37JC [1], where T is the temperature, V^s^ is the rate at which entropy 
is being generated in a unit volume and C is the bulk viscosity coefficient. The 
second law of thermodynamics can be written as TV^s'^ > => 3ffC ^ 0. 
Since H is positive in an expanding universe then C has to be positive in order to 



^We use this particular definition of S7m for computational conveniences. Pcritic 
3H^/8ttG is the critical density today, where Hq is the Hubble constant. 
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Figure 1: Joint confidence intervals for the parameters (CojHo) for a bulk vis- 
cous matter-dominated universe, spatially flat and without any dark energy 
component. We show the confidence intervals of 68.3%(lcr), 95.4%(3(t) and 
99.73%(5(t). The constraints are derived from the Gold 2006 SNe la data sam- 
ple alone (left), and the joint Gold 2006 -|- CMB -|- BAO cosmological tests 
(right). 

preserve the validity of the second law of thermodynamics. Thus, in the present 
model a necessary condition is Cm > 0. 

Constraining the bulk viscous cosmological model We compute the best 
estimated values and "the goodness-of-fit" of and Hq to the data through 
X^-minimization, using SNe alone, and also considering the CMB and BAO. 
Then we compute the confidence intervals for (^Qj ^^o) to constrain their possible 
values. We obtain as best estimates using only SNe la: (^o = 1-71, Ho = 
62.30) with a Xmin — 161.2 (Xd.o.f. = 0.895), where hereinafter Hq is in units 
of km • sec~i • Mpc~^ And using SNe la, CMB and BAO together: (Co = 
-0.30, ffo = 52.97) with a xLn = 400.5 (xi^.f. = 2.225). 



It can be seen that when we consider only SNe la we obtain Co ^ with a 
99.7% confidence level and a reasonable value for the Hubble constant in the 
interval 59.5 > Hq > 65.4 as well as for Xd o f = 0.895. However, when the same 
analysis using the three cosmological tests together (SNe la, CMB and BAO) is 
performed, we obtain negative values of Co with at least 99.7 % confidence level 
that disagree with the SLT, a not so good estimation for Hq of 51.4 > Hq > 54.6, 
and a bad xl o f ~ 2.225. These results are illustrated in Figure [1] 

Therefore, it is possible to conclude that this model is not a good candidate 
for explaining the observed accelerated expansion of the universe as suggested 
in [ill[5]. 
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